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Abstract 

Background: Low temperatures play key roles in the development of most plant diseases, mainly because of their 
influence on the expression of various virulence factors in phytopathogenic bacteria. Thus far, studies regarding this 
environmental parameter have focused on specific themes and little is known about phytopathogenic bacteria 
physiology under these conditions. To obtain a global view regarding phytopathogenic bacteria strategies in 
response to physiologically relevant temperature changes, we used DNA microarray technology to compare the 
gene expression profile of the model bacterial pathogen P. syringae pv. phaseolicola NPS3121 grown at 18°C 
and 28°C. 

Results: A total of 236 differentially regulated genes were identified, of which 133 were up-regulated and 103 were 
down-regulated at 18°C compared to 28°C. The majority of these genes are involved in pathogenicity and virulence 
processes. In general, the results of this study suggest that the expression profile obtained may be related to the 
fact that low temperatures induce oxidative stress in bacterial cells, which in turn influences the expression of iron 
metabolism genes. The expression also appears to be correlated with the profile expression obtained in genes 
related to motility, biofilm production, and the type III secretion system. 

Conclusions: From the data obtained in this study, we can begin to understand the strategies used by this 
phytopathogen during low temperature growth, which can occur in host interactions and disease development. 
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Background disease development [1,2]. Regarding environmental 
During their life cycles, phytopathogenic bacteria possess conditions, temperature is a key factor in most plant dis- 
an epiphyte growth stage during which they can grow eases, which are favored mainly by low temperatures 
and reproduce on the surface of a plant without causing [1,3,4]. The influence of low temperature on disease de- 
disease. However, when conditions are favorable, the velopment is predominantly due to the expression of 
bacteria enter a pathogenic stage that leads to disease various pathogenicity and/or virulence factors by the 
development. It is known that a complex interaction be- pathogens, which influences plant health. Several bacter- 
tween key factors must exist for the development of the ial phytopathogens, such as Pseudomonas syringae and 
disease in plants, represented as "disease triangle". This Erwinia sp., produce disease in their host plants in re- 
involves the interaction of a susceptible host, a virulent sponse to low temperature, which appears to be the cue 
pathogen, and environmental conditions favorable for for these phytopathogens to produce virulence factors, 

including toxins, cell-wall degrading enzymes, and ef- 
fector proteins [4]. Thus, low temperatures are an im- 
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One of the most common and economically important 
diseases is the bean disease (Phaseolus vulgaris L.) known 
as "halo blight" because it causes major field crop losses. 
This disease, caused by the bacterial pathogen P. syringae 
pv. phaseolicola, affects both the leaves and pods [5-7]. 
Cool temperatures (less than 25°C) favor disease develop- 
ment, a condition that also favors production of the major 
virulence factor of the pathogen, known as phaseolotoxin 
[8,9]. Phaseolotoxin is a non-host specific and chlorosis- 
inducing toxin that acts as a reversible inhibitor of the en- 
zyme ornithine carbamoyltransferase (OCTase; EC2. 1.3.3), 
which catalyzes the conversion of ornithine to citruline in 
the arginine biosynthetic pathway [10,11]. The production 
of phaseolotoxin by P. syringae pv. phaseolicola is regu- 
lated mainly by temperature and is optimally produced at 
18°C-20°C, whereas at 28°C (the optimal growth tem- 
perature for this bacterium), the toxin is not detected 
[8,9]. Genes whose products are involved in phaseolotoxin 
synthesis are clustered within of a chromosomal region 
known as the "Pht cluster", whose expression is also low 
temperature (18°C) dependent [12]. Thus, like other phy- 
topathogenic bacteria, low temperatures play an important 
role in P. syringae pv. phaseolicola for the development of 
halo blight disease. However, research related to this en- 
vironmental parameter has focused only on phaseolotoxin 
production and little is known about the bacterium physi- 
ology under these conditions. Likewise, studies performed 
in other phytopathogenic bacteria have focused on specific 
topics regarding low temperature function [4]. 

Global knowledge about the strategies used by these 
phytopathogens, in terms of temperature change which 
influences virulence stage and disease development, is 
very scarce and most of these studies have focused on 
animal pathogens where high temperature caused this 
effect [13,14]. Therefore, this study was undertaken with 
the objective to understand how phytopathogenic bac- 
teria, in particular the bacterial pathogen P. syringae pv. 



phaseolicola NPS3121, respond to temperature changes 
related to the development of the most of plant diseases. 

Results and discussion 

Low temperature (18°C) negatively affects the growth 
rate of P. syringae pv. phaseolicola NPS3121 

To obtain a global view regarding the strategies used by P. 
syringae pv. phaseolicola NPS3121 in response to physio- 
logically relevant temperature changes, we used DNA 
microarray technology. We compared gene expression 
profiles in the P. syringae pv. phaseolicola NPS3121 wild- 
type (wt) strain grown at 18°C and 28°C in M9 minimal 
media. These temperatures represent conditions that ei- 
ther favor the development of the disease (18°C) or do not 
(28°C) [8]. Initially, to evaluate the effect of temperature 
and establish the growth stage for this study, we per- 
formed bacterial growth curves of the P. syringae pv. 
phaseolicola NPS3121 strain grown under the conditions 
mentioned above. The results showed that at low tem- 
perature (18°C), the growth rate of the bacteria decreases 
approximately 0.5-fold relative to 28°C (Figure 1A). This 
behavior was reproducible in all performed kinetics. The 
effect of low temperature on the growth rate of several 
Pseudomonas syringae strains, including pv. phaseolicola, 
had been previously observed with similar results to this 
study [15]. Because previous results from our group indi- 
cated that during the transition phase, low temperature- 
induced differential expression in the phaseolotoxin syn- 
thesis genes (Pht cluster) occurs [12], we performed this 
study with cells harvested during this growth stage, which 
allowed us to use this cluster as a control for the micro- 
array and ensure the virulent stage of the bacterium. Thus, 
parallel cultures of P. syringae pv. phaseolicola NPS3121 
grown at 28°C and 18°C were harvested at the transition 
phase and RNA was extracted. The results presented in 
this work reflect the adapted state and significant genes, 
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whose expression is differentially maintained over long- 
term growth at a given temperature. 

Differentially expressed genes at 18°C are distributed 
throughout the chromosome and comprise several 
functional categories 

The differentially expressed genes were identified using 
a cut-off criteria of >1.5 for up-regulated and <0.6 for 
down-regulated genes (p-value < 0.05). A total of 236 dif- 
ferentially regulated genes were identified, of which 133 
were up-regulated and 103 were down- regulated at 18°C 
relative to 28°C. Analyses about the distribution and lo- 
cation of the genes in the P. syringae pv. phaseolicola 
1448A sequenced genome, showed that the differentially 
expressed genes at 18°C are not located in a single 
chromosomal region of P. syringae pv. phaseolicola, but 
rather are distributed throughout the genome. Further- 
more, only down-regulated genes were distributed in 
both plasmids of this strain (Figure 2). This pattern of 
distribution had been observed in preliminary assays, in 
which a Tn5-derived promoter probe was used to search 
for genes whose expression was temperature dependent; 
however, the authors reported the location of only a few 
genes throughout the genome [16]. 

For the purposes of this study, the differentially regu- 
lated genes were analyzed and manually grouped into 
categories based on their putative role in biological pro- 
cesses (Tables 1 and 2). In general, data analyses show 
that the majority of the differentially regulated genes 
relate to the pathogenicity and/or virulence process of 
the bacterium. 



Control and reliability of the microarray: Pht cluster 
genes and phaseolotoxin synthesis are induced at 18°C 

As expected, the microarray data showed that the Pht 
cluster genes (Cluster 1), involved in phaseolotoxin syn- 
thesis, were up-regulated at 18°C relative 28°C, thus 
confirming the previously reported data [12]. Thirteen of 
the 23 genes that comprise the Pht region were highly 
expressed at 18°C relative 28°C, which was consistent 
with the conditions of phaseolotoxin synthesis observed 
in the growth inhibition assays (Figure IB). Only 13 of 
the 23 Pht cluster genes were activated because only 
these genes are printed on our microarray. However, 
these genes represent the five transcriptional units that 
comprise the Pht region [12]. To validate the microarray 
data, one gene from each transcriptional unit was se- 
lected for validation of their expression pattern by 
RT-PCR analyses (Figure 3). The variability in Pht clus- 
ter gene expression levels observed could suggest diffe- 
rent regulation mechanisms for each of them. Thus far, 
it is known that there is transcriptional regulation for 
this group of genes mediated by temperature and only 
IHF protein has been identified as directly involved in 
the regulation of some of them [12,17]. The results re- 
garding the Pht cluster can also be used as control of 
the microarray, ensuring the reliability of the results ob- 
tained in this study. 

Genes involved in non-ribosomal peptide synthesis are 
induced at low temperature 

Another group of genes that was up-regulated at 18°C in 
P. syringae pv. phaseolicola NPS3121 comprise Cluster 2, 




Figure 2 Distribution and location of differentially expressed genes at 18°C in the P. syringae pv. phaseolicola genome. Differentially 
regulated genes were analyzed using the GenoMap software and their distribution and location in the bacterium genome was determined. The 
red bars depict the distribution of up-regulated genes and the green bars represent the down-regulated genes at 18°C. 



Arvizu-Gomez et al. BMC Microbiology 2013, 13:81 
http://www.biomedcentral.eom/1 471 -21 80/1 3/81 



Page 4 of 16 



Table 1 Genes up-regulated at 18°C in P. syringae pv. Table 1 Genes up-regulated at 18°C in P. syringae pv. 

phaseolicola NPS3121 phaseolicola NPS3121 (Continued) 



Gen/ORF 


Gene product 


Ratio 


PSPPH_ 


J 708 


ABC transporter, periplasmic 


1.84 


Cluster 1: Phaseolotoxin production (Pht cluster) 








substrate-binding protein 




PSPPH_ 


_4299 


Hypothetical protein {phtU) 


11.86 


PSPPH_ 


_2260 


UTP-glucose-1 -phosphate u ri dy ly It ra n sferase 


1.72 


PSPPH_ 


_4300 


Membrane protein, putative (phtT) 


8.70 


PSPPH_ 


_2542 


membrane protein, putative 


1.94 


PSPPH 


4301 


Adenylylsulfate kinase (phtS) 


1 3 ^n 

I JJU 


PSPPH_ 


_2643 


outer membrane efflux protein 


1.68 


PSPPH_ 


4302 


Conserved hypothetical protein (phtQ) 


D.Zj 


PSPPH 


2654 


lipoprotein, putative 


1.55 


PSPPH 


4305 


Hvnnthpfir^l nrnfpin fnhfO) 


8.78 


PSPPH_ 


_2842 


lipoprotein, putative 


1 .55 


PSPPH_ 


4306 


Hypothetical protein (phtM) 


15.90 


PSPPH_ 


_3226 


glycosyl transferase, group 1 family protein 


1 .72 


PSPPH 


4306 


Hypothetical protein (phtM) 


7.29 


PSPPH_ 


_3288 


predicted periplasmic lipoprotein 


2.05 


PSPPH 


4307 


n\/n i\/3fp nhn^nh^fp Hikin^p 

PEP/pyruvate binding subunit 


23.74 


PSPPH_ 
PSPPH_ 


_3810 
_3916 


lipoprotein 

membrane protein, putative 


1 .71 
1 .86 


PSPPH_ 


_4317 


Hypothetical protein 


11.52 


PSPPH_ 


_4139 


UDP-N-acetylglucosamine 


1.55 


PSPPH_ 


_4323 


Hypothetical protein 


2.13 






1 -carboxyvinyltransferase 




orgK 




control 


3.30 


PSPPH_ 


_4669 


acetyltransferase, GNAT family 


1.73 


pht A 




control 


4.96 


PSPPH_ 


_4682 


lipopolysaccharide biosynthesis 


3.03 


phtD 




control 








protein, putative 






6.50 














PSPPH_ 


_5220 


inner membrane protein, 60 kDa 


1.61 


desl 




control 


14.97 
















Cluster 5: Genes involved in motilitv 




phtL 




control 


7.64 












control 




PSPPH_ 


_0730 


type IV pi I us-associated protein, putative 


1.74 


phtMN 




1.81 












control 




PSPPH_ 


_0818 


type IV pilus prepilin peptidase PilD 


1.53 


omtA 




10.34 
















PSPPH_ 


_0820 


type IV pilus biogenesis protein PilB 


1.53 


Cluster 2: Genes involved in Non-ribosomal synthesis 










PSPPH_ 




transposon Tn7-like transposase protein A 




PSPPH_ 


J 200 


pili assembly chaperone 


2.03 


_4538 


1.67 








PSPPH_ 




transposon Tn7-like transposase protein B 




PSPPH_ 


_3387 


flagellar regulator FleQ 


1.64 


_4539 


1.70 












hypothetical protein PSPPH_4544 




PSPPH_ 


_3880 


CheW domain protein WspB 


1.94 


PSPPH_ 


_4544 


8.08 








PSPPH_ 






PSPPH_ 


_3881 


methyl-accepting chemotaxis protein WspA 


1.5 


_4544 


hypothetical protein PSPPH_4544 


7.26 








PSPPH_ 






Cluster 6: Oxidative stress response qenes and iron metabolism 


_4546 


hypothetical protein 


9.44 
















PSPPH_ 


J 309 


cysteine desulfurase IscS 


1.93 


PSPPH_ 


_4549 


hypothetical protein PSPPH_4549 


9.13 












major facilitator family protein 




PSPPH_ 


_1 31 1 


iron-sulfur cluster assembly protein IscA 


1.69 


PSPPH_ 


_4553 


15.83 








PSPPH_ 








PSPPH_ 


J 909 


RNA polymerase sigma-70 


1.59 


_4554 


arginine aminomutase, putative 


12.67 






family protein. pvdS 




PSPPH_ 


_4555 


conserved hypothetical protein 


/AO 


PSPPH_ 


J 923 


pyoverdine sidechain peptide 


1.70 


Cluster 3: Type VI secretion system 








synthetase I, epsilon-Lys module 




PSPPH_ 


_0122 


hep 


2.13 


PSPPH_ 


_2117 


FecR protein superfamily 


2.05 


PSPPH_ 


_0124 


hypothetical protein 


1.66 


PSPPH_ 


_3007 


iron ABC transporter, permease 


1.61 












protein, putative 




PSPPH_ 


_0125 


icmF 


1.94 


PSPPH_ 


_3274 


catalase KatB 


2.31 


PSPPH_ 


_0131 


HsiG 


1.61 


PSPPH_ 


_3274 


catalase KatB 


1.65 


PSPPH_ 


_0135 


hypothetical protein 


1.64 
















PSPPH_ 


_3753 


siderophore biosynthesis protein 


1.57 


PSPPH_ 


_4978 


prophage PSPPH06, putative reverse 


1.65 












transcriptase/maturase 




Cluster 7: Unknown function 




PSPPH_ 


_4979 


prophage PSPPH06, putative reverse 


2.68 


PSPPH_ 


_0317 


conserved hypothetical protein 


1.81 






transcriptase/maturase 




PSPPH_ 


_0611 


conserved hypothetical protein 


2.15 


PSPPH_ 


_4984 


prophage PSPPH06, site-specific 


1.70 


PSPPH_ 


_0612 


hypothetical protein PSPPH_0612 


1.94 






recombinase, phage integrase family 


















PSPPH_ 


J 142 


hypothetical protein PSPPHJ142 


1.58 


Cluster 4: Genes involved in membrane synthesis 














leucine-rich repeat domain protein 




PSPPH_ 


J 230 


hypothetical protein PSPPHJ230 


1.59 


PSPPH_ 


_1430 


1.75 








PSPPH_ 






PSPPH_ 


J 243 


conserved hypothetical protein 
hypothetical protein 


1.98 


_1464 


lipoprotein, putative 


1.66 


PSPPH_ 


J 637 


2.25 
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Table 1 Genes up-regulated at 18°C in P. syringae pv. 



phaseolicola NPS3121 (Continued) 

PSPPHJ835 conserved hypothetical protein 1.72 

PSPPHJ938 conserved hypothetical protein 1.52 

PSPPH_2103 conserved hypothetical protein 1.89 

PSPPH_2116 conserved hypothetical protein 1.58 

PSPPH_2147 hypothetical protein PSPPH_2147 2.51 

PSPPH_2578 conserved hypothetical protein 1.53 

PSPPH_3212 conserved hypothetical protein 1.53 

PSPPH_3262 conserved hypothetical protein 1.60 

PSPPH_5014 conserved hypothetical protein 1.52 
Cluster 8: Uncharacterized Function 

PSPPH_021 0 DNA repair protein RadC 1 .56 

PSPPH_0398 glutamate synthase, large subunit 2.63 

PSPPH_0581 radical SAM domain protein 1.53 

PSPPH_0620 DNA primase 2.48 

PSPPH_0622 O-sialoglycoprotein endopeptidase 1.87 

PSPPH_0625 2-amino-4-hydroxy-6- 1 .62 
hydroxymethyldihydropteridine 
pyrophosphokinase 

PSPPH_0627 SpoVR like family protein 2.1 0 

PSPPH_0629 protein kinase 1.65 

PSPPH_0703 phosphonate ABC transporter 1 .73 
permease protein phnE 

PSPPHJ 1 41 ISPsy20, transposase IstB 1 .5 1 

PSPPH_1 1 50 conserved domain protein-Divergente 1.61 
HU family 

PSPPHJ 1 79 DNA-binding response regulator GltR 1 .54 

PSPPHJ 244 transcriptional regulator, AsnC family 1.89 

PSPPHJ 306 RNA methyltransferase, TrmH family, group 1 1 .545 

PSPPHJ 378 Methionyl-tRNA synthetase 2.56 
(Methionine-tRNA ligase)(MetRS) 

PSPPHJ 406 ATP-dependent helicase, DinG family 1.77 

PSPPHJ 468 nucleic acid binding protein 1.58 

PSPPHJ 595 transcriptional regulator, GntR family 2.58 

PSPPHJ 661 cvpA family protein 1 .66 

PSPPHJ 746 oxidoreductase, aldo/keto reductase family 1 .92 

PSPPH_2216 zinc carboxypeptidase domain protein 1.89 

PSPPH_2221 precorrin-4 C1 1 -methyltransferase 1.52 

PSPPH_2506 L-arabinose ABC transporter, periplasmic 1.62 
L-arabinose-binding protein 

PSPPH_2551 oxidoreductase, putative 1.84 

PSPPH_2563 transcriptional regulator, GntR family 1.53 

PSPPH_2580 transcriptional regulator, LysR family 1.97 

PSPPH_2620 5-methyltetrahydrofolate-homocysteine 1.85 
methyltransferase 

PSPPH_2690 oxidoreductase, FAD-binding, putative 1.56 

PSPPH_2781 TspO/MBR family protein 1 .99 

PSPPH_2840 sodium/hydrogen exchanger family protein 1.55 



Table 1 Genes up-regulated at 18°C in P. syringae pv. 
phaseolicola NPS3121 (Continued) 



PSPPH 


2847 


npnpr^l ^prrptinn nathwav 
protein GspK, putative 


1 .89 


PSPPH 


_3045 


transporter, AcrB/AcrD/AcrF family 


1.64 


PSPPH 


_3252 


glycolate oxidase, GlcD subunit 


1.97 


PSPPH 


_3291 


oxidoreductase, molybdopterin-binding 


1.88 


PSPPH 


_3294 


DNA-binding heavy metal response regulator 


1.81 


PSPPH 


_3654 


transcriptional regulator, TetR family 


1.51 


PSPPH 


_3906 


sensor histidine kinase 


1.65 


PSPPH 


_3946 


DNA repair protein RecO 


1.66 


PSPPH 


_3962 


DNA-binding response regulator TctD 


1.77 


PSPPH 


_4137 


histidinol dehydrogenase 


1.63 


PSPPH 


_4151 


RNA polymerase sigma-54 factor RpoN 


1.69 


PSPPH 


41 52 


llUUoUllldl oUUUIIIL 1 1 1 Lcl I dLc piULclll 


1 86 


PSPPH. 


_4332 


DNA repair protein RadA 


1.76 


PSPPH 


_4372 


RNA 2'-phosphotransferase 


1.55 


PSPPH 


_4634 


bmp family protein 


2.99 


PSPPH 


_4641 


YccA 


1.68 


PSPPH 


_4717 


dethiobiotin synthetase 


2.09 


PSPPH 


_4866 


proline-specific permease proY 


1.54 


PSPPH 


_4925 


imidazole glycerol phosphate synthase, 
glutamine amidotransferase subunit 


1.62 


PSPPH 


_5142 


oxaloacetate decarboxylase alpha subunit 


2.35 



The described functions were obtained from the literature. The up-regulated 
genes were identified using cutoff criteria >1.5 of ratio. The ratio is in relation 
to expression levels obtained between 18°C and 28°C (18°C/28°C). Control: 
corresponds to genes obtained by PCR amplification that were printed in 
the microarray. 



corresponding to genes involved in non-ribosomal peptide 
synthesis (NRPS) [18]. NRPS is an alternative pathway that 
allows production of polypeptides via a different mecha- 
nism than the traditional translation pathway. Peptides 
are created by enzymatic complexes called synthetases. 
Through NRPS, some bacteria produce several secondary 
metabolites, such as siderophores, antibiotics, or toxins 
that contribute to the fitness and/or pathogenicity of the 
bacterium [19,20]. In our microarray, six genes that en- 
code four hypothetical proteins (PSPPH_4544, PSPPH_ 
4546, PSPPH_4549, and PSPPH_4555), a facilitator family 
protein (PSPPH_4553), and an arginine aminomutase 
(PSPPH_4554) were highly up-regulated at 18°C. These 
genes are located in a 27 kbp fragment, which also en- 
codes a polyketide synthetase domain protein (PSPPH_ 
4547) and a non-ribosomal peptide synthetase (PSPPH_ 
4550). This region is delimited by genes encoding for 
transposases (PSPPH_4538 and PSPPH_4559), which were 
also induced in our microarray (Table 1, Cluster 2). Of all 
the genes of this region, only six genes were printed on 
the microarray and all of these were induced in the condi- 
tions evaluated. RT-PCR analyses confirmed the profile 
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Table 2 Genes down-regulated at 18°C in P. syringae pv. 
phaseolicola NPS3121 

Gen/ORF Gene product Ratio 



Table 2 Genes down-regulated at 18°C in P. syringae pv. 



phaseolicola NPS3121 (Continued) 

Cluster 13: Genes related with nucleic acids synthesis 

PSPPH_4598 DNA-directed RNA polymerase, beta' subunit 0.59 

PSPPH_4599 DNA-directed RNA polymerase, beta' subunit 0.57 

PSPPH_2495 DNA polymerase II 0.57 

PSPPH_B0043 DNA topoisomerase III 0.64 

PSPPH_A0002 Replication protein 0.54 
Cluster 14: Unknown function 

PSPPH_0220 conserved hypothetical protein 0.64 

PSPPH_0609 hypothetical protein PSPPH_0609 0.54 

PSPPH_2482 conserved hypothetical protein 0.63 

PSPPH_2855 hypothetical protein PSPPH_2855 0.43 

PSPPH_3333 conserved hypothetical protein 0.36 

PSPPH_3625 conserved hypothetical protein 0.59 

PSPPH_4047 conserved hypothetical protein 0.66 

PSPPH_A0040 hypothetical protein PSPPH_A0040 0.66 

PSPPH_B0048 conserved hypothetical protein 0.60 
Cluster 15: Uncharacterized function 

PSPPH_001 2 glycyl-tRNA synthetase, alpha subunit 0.63 

PSPPH_0033 3-oxoadipate enol-lactonase, putative 0.65 

PSPPH_0072 membrane protein, putative 0.63 

PSPPH_0080 ATP-dependent DNA helicase Rep 0.43 

PSPPH_01 1 7 phospholipase D family protein 0.63 

PSPPH_0215 aldehyde dehydrogenase family protein 0.35 

PSPPH_0296 colicin/pyocin immunity family protein 0.58 

PSPPH_0360 periplasmic glucan biosynthesis protein 0.63 

PSPPHJ 072 oxidoreductase, short chain 0.57 
dehydrogenase/reductase family 

PSPPHJ 181 glucose ABC transporter, periplasmic 0.65 
glucose-binding protein, putative 

PSPPH_1 21 1 cytochrome o ubiquinol oxidase, subunit I 0.55 

PSPPHJ 508 acetyltransferase, GNAT family 0.35 

PSPPHJ 5 1 8 ATP-dependent DNA helicase RecQ 0.53 

PSPPHJ 575 CAIB/BAIF family protein 0.65 

PSPPHJ 759 plasmid stabilization system family protein 0.53 

PSPPHJ 762 transcriptional regulator, AsnC family 0.54 

PSPPHJ 91 7 cation ABC transporter, periplasmic 0.60 
cation-binding protein 

PSPPHJ 921 peptidase 0.58 

PSPPHJ 963 electron transfer flavoprotein-ubiquinone 0.38 
oxidoreductase, putative 

PSPPH_2053 membrane protein, putative 0.65 

PSPPH_2057 2-methylcitrate synthase 0.62 

PSPPH_2159 dehydrogenase, isocitrate/isopropylmalate family 0.60 

PSPPH_2246 4-alpha-glucanotransferase 0.66 

PSPPH_2695 peptide ABC transporter, 0.45 
permease protein 



Cluster 9: Alginate synthesis 



PSPPHJ 112 


alginate biosynthesis protein AlgX 


0.52 


PSPPHJ 113 


alginate biosynthesis protein AlgG 


0.19 


PSPPHJ 114 


alginate biosynthesis protein AlgE 


0.18 


PSPPHJ 115 


alginate biosynthesis protein AlgK 


0.19 


PSPPHJ 118 


alginate biosynthesis protein AlgD 


0.46 


PSPPHJ 119 


conserved hypothetical protein 


0.46 


olgD 


olgD (control) 


0.25 


Cluster 10: Plant-Pathogen interactions 




PSPPH_A0075 


type III effector HopW1-2, truncated 


0.60 


PSPPH_A0127 


type III effector HopABI 


0.42 


PSPPH_A0127 


type III effector HopABI 


0.65 


PSPPH_A0127 


virA type III HopABI (control) 


0.57 


PSPPH_A0120 


avrC type III effector AvrB2 (control) 


0.53 


PSPPH_A0010 


avrD type III effector hopD1 (control) 


0.56 


PSPPH_3992 


pectin lyase 


0.62 


PSPPH_3993 


acetyltransferase, GNAT family 


n ^7 

U.D / 


PSPPH_A0072 


polygalacturonase 


0.50 


Cluster 1 1 : Type IV secretion system 




PSPPH_B0022 


transcriptional regulator, PbsX family 


0.65 


PSPPH_ 
B0023 


transcriptional regulator 


0.64 


PSPPH_ 

B0025 


conjugal transfer protein 


0.65 


PSPPH_ 

B0027 


conjugal transfer protein 


0.65 


PSPPH_ 

B0028 


conjugal transfer protein 


0.61 


PSPPH_ 
B0031 


conjugal transfer protein 


0.65 


PSPPH_ 
B0032 


conjugal transfer protein 


0.61 


PSPPH_ 
B0034 


conjugal transfer protein 


0 62 


PSPPH_ 
B0035 


conjugal transfer protein 


0.66 


PSPPH_ 

B0036 


conjugal transfer protein 


0.51 


PSPPH_ 
B0041 


conjugal transfer protein 


0.58 


Cluster 12: Heat-shock proteins 




PSPPH_0381 


heat shock protein HsIVU, ATPase subunit HsIU 


0.65 


PSPPH_0742 


dpB protein 


0.54 


PSPPH_4077 


chaperonin, 60 kDa. groEL 


0.29 


PSPPH_4206 


dnaK protein 


0.28 


PSPPH_4206 


dnaK protein 


0.57 


PSPPH_4207 


heat shock protein GrpE 


0.65 
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Table 2 Genes down-regulated at 18°C in P. syringae pv. 
phaseolicola NPS3121 (Continued) 



PSPPH 


2868 


m^inr farilit^tnr familv tr^n^nnrfpr 
1 1 icijvji i ul.ii i La laiiiiiy li cii lopLJi lci 


0.63 


PSPPH 


2892 


TnnR-HpnpnHpnt ^iHprnnhnrp rprpntnr ni it^tk/p 


0.62 


PSPPH 


2897 


x/PKini^h^rtin nnn-rihn^nm^l npntiHp ^vnthpf^p 
yci on iiakja^-Lii i i i i ikj^o^i i lai pcpuuc oyi ili iclcoc 


0.40 


PSPPH 


2899 


vprQinishsrtin nnl\/kptiHp/nnn-ri hn^nins 1 

ycrijii iiakja^-Lii i uj^jiy r\cr li cy i iui i iikj^jowiiiai 

peptide synthetase 


0.58 


PSPPH_ 


_2904 


isochorismate synthase 


0.55 


PSPPH_ 


_3100 


isocitrate dehydrogenase, NADP-dependent 


0.63 


PSPPH. 


_3251 


maleylacetoacetate isomerase 


0.53 


PSPPH. 


_3528 


acetate-CoA ligase 


0.52 


PSPPH. 


_3558 


aconitate hydratase 2 


0.61 


PSPPH. 


_3782 


porin D 


0.42 


PSPPH. 


_3985 


3-oxoacyl-[acyl-carrier protein] reductase 


0.54 


PSPPH. 


_4221 


unnamed protein product 


0.44 


PSPPH. 


_4654 


smtA protein 


0.47 


PSPPH. 


_4703 


coenzyme PQQ biosynthesis protein PqqF 


0.32 


PSPPH. 


_4805 


oxidoreductase FAD-binding domain/ 
oxidoreductase NAD-binding domain/2Fe-2S 

iron-^i ilfi ir Hi Ktpr hinHinn domain nrntpin 

mui i juii ui L-iuoLcri kj\ i i \A 1 1 ly ^ ui i ia 1 1 i ij\ \j lci i i 


0.55 


PSPPH 


4833 


Rhs family protein 


0.33 


PSPPH 


4859 


tran^nnrtpr Rf~f~T familv 

Liaiio^jwiLLri, i i a i i 1 1 ly 


0.65 


PSPPH 


4869 


cadmium-translocating P-type ATPase 


0.54 


PSPPH 


4885 


D-^-nhn^nhnnK/rpr^tp r\p)r\\/r\rrir\pr\?\<^p 
l> j pi iLjopi iuy ly L.CI a lc uci lyuiLjyci iojc 


0.56 


PSPPH 


4938 


^minn ^riH ARf~ Ir^n^nnrtpr ATP-hinHinn nrntpin 

Ul 1 III ILJ UL_IU /\L>^_ LI Ql loLJLJI LCI , nl T kJIIIUIIiy LJI LJLCII 1 


0.61 


PSPPH. 


_4962 


prophage PSPPH06, C4-type zinc finger protein, 
DksA/TraR family 


0.35 


PSPPH 


5024 


acetyltransferase, GNAT family 


0.64 


PSPPH 


5027 


ar pt\/|trancfpracp (^.KIAT f^milv 
aL.CLyiLiaiioicia.jc, ljin/\i laiiiiiy 


0.64 


PSPPH 


51 70 


acy 1 LI al lol cl aoc lalllliy |JIULclll 


0 60 


PSPPH 


A0062 


1 v^R-fpi mi l\/ tr^ rKrri ntinn rpniil^tnr ^inR 
i_yji\ laiiiiiy li ai loci iljlilji i icyuiaLUi jii in 


0.45 


PSPPH. 


_A0083 


IS801, transposase 


0.64 


PSPPH. 


_A0109 


sulfotransferase, putative 


0.49 


PSPPH. 


_A0129 


Yersinia/Haemophilus virulence 
surface antigen family 


0.53 


PSPPH. 


_A0132 


ISPsyl 6, transposase 


0.66 


PSPPH. 


_A0145 


conjugal transfer protein 


0.56 


PSPPH. 


_B0004 


RulB protein 


0.63 


PSPPH. 


_B0050 


relaxase, putative 


0.65 


PSPPH. 


_B0059 


exeA-like protein 


0.64 



The described functions were obtained from the literature. The down-regulated 
genes were identified using cutoff criteria < 0.6 of ratio. The ratio is in relation 
to the expression levels obtained between 18°C and 28°C (18°C/28°C). Control: 
corresponds to genes obtained by PCR amplification that were printed in 
the microarray. 



expression and temperature dependence of PSPPH_4553 
and PSPPH_4554, thus validating the microarrays results 
(Figure 3). The expression of at least one of these genes 
(PSPPH_4550) in temperature dependence had been pre- 
viously observed with similar results [21]. In P. syringae 







RT-PCR Expression 




18°C 28°C 


Phaseolotoxin synthesis 




argK 




phtA 




desl 


MB^K 


phtL 




amtA 




Non-ribosomal peptide synthesis 




Hypothetical protein PSPPH_4547 




Hypothetical protein PSPPH_4553 




Hypothetical protein PSPPH_4554 




Motility 




FleQ 




Oxidative stress and iron metabolism 




r\iNM [juiyi i ici doc s>iyiiia-/u idiimy piuieiii pvuo 


1 6 


Catalase KatB 


1W 2.3 


Plant-pathoQen Interaction 




Type 1 1 1 effector virPphA 


0 5 


Type III effector avrPphC 


____ n r 


Type III effector avrPphD 




Pectin lyase pnIA 




Polygaracturonase 


0.5 


Alginate synthesis 




algD 




Uncharacterized function 




RNA Polymerase sigma-54 factor rpoN 




Type VI Secretion System 




PSPPH_0122 hep 




PSPPH_0127 hsiJ 




PSPPH_0131 hsiG 


■■■■■ 


16S ribosomal 




Figure 3 Microarray validation using RT-PCR analyses. RT-PCR 


validates the microarray results. a Corresponds to expression levels 


obtained in the microarray for these genes. The remaining genes do 


not show expression levels because they were not printed on 


the microarray. 





pv. phaseolicola NPS3121, it has been suggested that 
NRPS genes are part of a genomic island (IG) acquired by 
horizontal transfer and is postulated to be involved in 
phaseolotoxin synthesis during peptide assembly. How- 
ever, only the PSPPH_4550 gene has been demonstrated 
to have a role in this process [21], Based on this hypoth- 
esis, the profile expression obtained for this group of 
genes at 18°C could be congruent with the differential ex- 
pression of the Pht cluster genes and the conditions for 
phaseolotoxin synthesis. However, the RT-PCR results for 
the PSPPH_4547 gene showed that the expression of this 
gene is independent of temperature, presenting constitu- 
tive behavior at both temperatures (Figure 3). Knowledge 
regarding the role of this P. syringae pv. phaseolicola gene 
group is limited and experimental work is still necessary. 
Likewise, is necessary to evaluate whether there is a rela- 
tionship between these genes and phaseolotoxin synthesis 
genes, as has been previously proposed, or whether these 
genes participate in different biological processes that con- 
tribute to the fitness of the bacterium in low temperatures. 
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In P. syringae pv. phaseolicola NPS3121, the Type VI 
secretion system (T6SS) is regulated by temperature 

Recently, a new secretion system has been recognized, 
called the Type VI secretion system (T6SS). This system 
is encoded within the genomes of most Gram negative 
bacteria, including plant, animal, and human pathogens, 
as well as environmental strains. The T6SS components 
are usually encoded by a gene cluster that is thought to 
form a genomic island whose composition and number 
varies among species [22-24]. The in silico analyses have 
revealed that the genome of P. syringe pv. phaseolicola 
1448A carries only one putative T6SS gene cluster (HSI) 
that comprises the region from PSPPH_0119 to PSPPH_ 
0135. Furthermore, several genes putatively encoding 
some proteins of this system are scattered in the genome 
of this bacterium [24]. The microarrays results showed 
the induction of eight genes encoding proteins putatively 
involved in the T6SS in P. syringe pv. phaseolicola 
NPS3121 (Cluster 3). The PSPPH_0122 gene encodes a 
hemolysin-coregulated (Hep) protein homolog, in addition 
to be an essential component of the secretion machinery, 
acts as an effector protein that is secreted through this 
system. The PSPPH_0124 gene encodes a hypothetical 
protein and the PSPPH_0125 gene encodes the IcmF 
protein, which in conjunction with the DotU protein 
(PSPPH_0126), act as associated structural proteins that 
anchor the secretion system in the cell membrane [25]. 
Within this cluster is also the PSPPH_0131 gene encoding 
the hsiG protein and the PSPPH_0135 gene that encodes 
a hypothetical protein. The expression of these genes 
at low temperature was confirmed by RT-PCR (Figure 3). 
Although PSPPH_ 4978, PSPPH_ 4979, and PSPPH_ 
4984, which encode prophage PSPPH06 proteins, are not 
involved in T6SS, these genes were include within this 
group because their adjacent genes (PSPPH_ 4980 and 
PSPPH_ 4985) putatively encode Hep proteins [24], which 
may be responsible for the induction levels obtained. This 
finding is being evaluated in our laboratory. The T6SS has 
been shown to play a key role in the virulence and patho- 
genesis of diverse bacterial pathogens, in some cases, by 
the secretion of effector proteins or toxins. However, its 
complete mechanism of action is poorly understood. The 
function of this system is not restricted to pathogenic pro- 
cesses because the T6SS also participates in other pro- 
cesses such as biofilm formation, stress sensing, symbiosis, 
root colonization, and nodule formation [26,27]. The role 
of the putative T6SS gene cluster in P. syringae pv. 
phaseolicola NPS3121 has not been evaluated so more ex- 
perimental work is required. However, it has been demon- 
strated that T6SS in P. syringae pv. syringae B728a, which 
is phylogenetically identical to P. syringae pv. phaseolicola 
T6SS, it is not essential for leaf colonization and develop- 
ment of the disease [28]. Several reports have demon- 
strated that expression of the T6SS gene cluster is tightly 



regulated in different environmental conditions and low 
temperatures contribute to the expression of these genes 
in some pathogens [29]. This phenomenon is similar to 
our observation that low temperature (18°C) regulates 
T6SS genes expression. To our knowledge, this is the first 
report about expression of these genes of P. syringae pv. 
phaseolicola NPS3121 and the influence of low temperature 
on their expression. 

Cell envelope-associated changes are induced by low 
temperature 

A universal response to low temperature includes changes 
in the lipid composition of membranes to help cope with 
the decrease in membrane fluidity caused by the cold. 
Microorganisms respond by increasing the unsaturated 
fatty acids level in membrane phospholipids, which helps 
to maintain membrane homeoviscosity so that its function 
is not affected. There are a variety of mechanisms that can 
alter membrane phospholipid composition in response to 
temperature change [30]. The conversion of saturated 
fatty acids into unsaturated fatty acids by desaturases en- 
zymes is one of these pathways [30,31]. In our microarray 
and RT-PCR analyses (Figure 3, Cluster 1), the desl gene 
encoding a fatty acid desaturase was induced at 18°C, 
which might be involved in the unsaturation process, in a 
similar manner to the reported desA and des genes from 
Synechosysteis sp. PCC6803 and Bacillus subtilis, respec- 
tively. It has been observed that deletion of the des gene in 
B. subtilis produces a cold-sensitive phenotype and slower 
growth, thus demonstrating its role during adaptation to 
low temperatures [32]. In P. syringae pv. phaseolicola 
NPS3121, the desl gene is part of the phtD transcriptional 
unit of the Pht cluster (Cluster 1). Mutation of this gene 
produces a non-toxigenic phenotype relative to the wt 
strain. However, the relationship of desl with phaseolotoxin 
synthesis is still unknown [12]. Additionally, it has been 
observed that mutation in the desl gene decreases the 
growth rate at 18°C relative to the wt strain, suggesting a 
cold-sensitivity in the mutant strain (unpublished data). 

Another of the mechanisms reported to be involved in 
membrane lipid composition changes correspond to de 
novo synthesis. The fabF and IpxP genes induced by low 
temperature participate in this process [33]. (3-ketoacyl- 
ACP synthase II, the fabF gene product, converts 
palmitoleic acid to ds-vaccenic acid, which is in turn 
transferred by an acyltransferase (LpxP) into lipid A, a 
component of polysaccharides [33,34]. Although these 
two genes were not found in our microarray, several 
genes involved in cell wall biogenesis and membrane 
synthesis were identified (Cluster 4). These include the 
murA gene (PSPPH_4139) that is involved in peptidogly- 
can synthesis (a major component of cell wall), the 
PSPPH_4682 gene involved in polysaccharide synthesis, 
as well as three genes PSPPH_4669, PSPPH_3226, and 
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galU (PSPPH_2260) that encode an acetyl-, glycosyl- 
and uridyl- transferase, respectively, which are likely 
associated with the transfer of these groups during poly- 
saccharides synthesis. Additionally, it has been demon- 
strated that during cell envelope biogenesis, there is an 
increase in outer membrane lipoproteins, which increase 
connections with the cell wall [34,35]. In our analyses 
four genes (PSPPH_ 1464, PSPPH_2654, PSPPH_2842, 
and PSPPH_3810) encoding lipoproteins were induced, 
which may be related to outer membrane synthesis. The 
microarray results suggest that membrane component 
synthesis is activated in the conditions of our study and 
these changes are likely related to cell envelope remodel- 
ing to adapt to low temperatures. 

Low temperature induces expression of motility genes in 
P. syringae pv. phaseolicola NPS3121 

Cluster 5 comprises genes induced at 18°C that are in- 
volved in bacterium motility. The data suggest that 
chemotaxis and rotation of flagella processes function 
in low temperatures on P. syringae pv. phaseolicola 
NPS3121. Two genes, PSPPH_3880 that encodes the 
membrane-bound methyl accepting chemotaxis protein 
(MCP)-like receptor WspA, and PSPPH_3881, that en- 
codes the CheW-like scaffolding protein WspB, showed 
high transcripts levels at 18°C relative 28°C (Table 1). 
WspA and WspB are related to the chemotaxis process. 
Chemotaxis, as well as other types of taxis (e.g., thermo- 
taxis), enables bacteria to approach beneficial environ- 
ments and escape from hostile ones. Depending on the 
parameter monitored, bacteria will respond by either 
swimming toward attractants or retreating from repellants. 
Thus, the signal sensed by chemotaxis causes changes in 
flagellum motility [36]. Flagella gene expression is acti- 
vated by the FleQ protein, which appears to be the highest 
regulator in the hierarchical flagellar biogenesis regulatory 
cascade in P. aeruginosa [37]. FleQ (PSPPH_3387) was in- 
duced in our study at 18°C and its expression was vali- 
dated by RT-PCR (Figure 3), suggesting that the motility 
of P. syringae pv. phaseolicola NPS3121 is favored under 
this condition. Furthermore, four genes related to pili for- 
mation, which is also involved in bacterial movement, were 
induced at low temperature: PSPPH_0730 that encodes 
type IV pilus-associated protein, PSPPH_1200 that en- 
codes a pili assembly chaperone, PSPPH_0818 that en- 
codes PilD protein, and PSPPH_0820 that encodes PilB 
protein. Each of these genes has been associated with 
P. syringae pv. phaseolicola virulence because of their role 
in adhesion to the surface of host plants to initiate infec- 
tion [38] . It has been reported that RpoN sigma factor reg- 
ulates the expression of genes required for pili and flagella 
biosynthesis in P. aeruginosa [37,39]. Our microarray data 
and RT-PCR assays showed that the PSPPH_4151 gene 
(Cluster 8), which encodes the RpoN protein, was induced 



at 18°C, suggesting a similar regulation may occur in our 
strain (Figure 3). The results obtained suggest that 
P. syringae pv. phaseolicola NPS3121 motility is regulated 
by temperature, similar to those observed in the pathogens 
Helicobacter pylori and E. colU whose motility patterns are 
altered by temperature changes [33,40]. 

To assess whether these changes in the gene expres- 
sion generate a motility phenotype in P. syringae pv. 
phaseolicola related to temperature, we evaluated the 
motility pattern of this bacterium at 18°C and 28°C. To 
ensure that the bacteria were in the same physiological 
condition as when the microarray analysis was performed, 
P. syringae pv. phaseolicola NPS3121 cells grown at 18°C 
and 28°C (OD 600 : 1.1 and 1.2) were inoculated in semi- 
solid M9 media containing 0.3%, 0.4%, and 0.5% agar and 
incubated at 18°C and 28°C. The results showed that 
under these conditions the bacterium was not motile des- 
pite gene induction at 18°C (Figure 4A). Additionally, mo- 
tility assays in KB media were performed, using the 
conditions that have demonstrated motility in related 
pathovars [41,42]. Plates and glass tubes with semisolid 
KB media were used to evaluate motility at the mentioned 
temperatures. Again, the P. syringae pv. phaseolicola 
NPS3121 strain was not motile under these conditions 
compared to P. syringae pv. tomato DC3000 and P. 
syringae pv. tabaci, which showed motility at both temper- 
atures and where it was observed that low temperatures 
appear to affect their motility (Figures 4B and 4C). This 
non-motile phenotype of P. syringae pv. phaseolicola 
NPS3121 had been previously reported [41,43], and fur- 
ther experiments are required to determine the conditions 
in which this bacterium can be motile and to evaluate the 
effect of low temperature in this process. 

Low temperature induces oxidative stress and iron 
metabolism 

Another group of genes differentially expressed at 18°C 
correspond to genes related to iron metabolism (Cluster 6). 
Iron fulfills a vital role in virtually all organisms because 
of its participation in several cellular processes. Because 
iron is in short supply in many habitats, bacteria secrete 
siderophores, compounds that are specific iron chelators, 
to mobilize inside the cell through membrane receptor 
molecules [44]. Two genes, PSPPH_3753 that encodes a 
protein related to siderophore synthesis and PSPPH_1923 
that is involved in pyoverdine synthesis (a major sidero- 
phore of the fluorescent Pseudomonas sp), were induced at 
18°C relative to 28°C [45]. Likewise, the gene encoding 
sigma factor protein PvdS, which is required for expression 
of pyoverdine synthesis genes, was induced under these 
conditions [46]. The induction of this PvdS protein was 
validated by RT-PCR analysis (Figure 3). One gene encod- 
ing the regulatory protein FecR (PSPPH_2117) and proteins 
involved in iron transport were also included in this group. 
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NPS3121 DC3000 



Figure 4 Motility Tests of the P. syringae pv. phaseolicola NPS3121 strain. (A) Motility assays performed in semisolid M9 media containing 
0.3%, 0.4%, and 0.5% agar at 18°C and 28°C. (B) Motility assays in semisolid KB media containing 0.3% (left) and 0.5% (right) agar. (C) The results 
obtained using the stab technique in M9 and KB media. 



It is known that in P. aeruginosa, the Fur protein is 
the master regulator of iron homeostasis. It represses 
pyoverdine synthesis via negative regulation of the pvdS 
gene under high iron concentrations. However, in iron- 
limiting conditions, Fur repression is released and tran- 
scription can occur [47]. It has been reported that PvdS 
sigmulon is conserved among the fluorescent pseudomo- 
nads, including the P. syringae group [46]. Although the fur 
gene was not printed on our microarray, the functional 
status of Fur protein can be inferred as inactive because 
the genes regulated by this protein are induced in the 



conditions evaluated. This expression profile simulates con- 
ditions of iron deficiency. 

To phenotypically evaluate changes in the expression of 
siderophores synthesis genes in function of temperature, 
we performed quantitative analyses of siderophores at 
18°C and 28°C. The results of these assays showed that 
at 18°C, the amount of siderophores in the culture su- 
pernatant was higher (58.6 ± 0.39 \iM) compared to when 
the bacterium is grown at 28°C (20.53 ± 0.844 \sM). Thus, 
the results demonstrate that low temperature induces 
siderophores production by the bacterium. 
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Additionally, it has been reported that in several bac- 
teria, the Fur protein positively regulates the expression 
of genes involved in various pathways in response to 
large iron amounts, such as oxidative stress genes (e.g. 
catalases) [47]. In our microarray, the PSPPH_3274 gene 
(encoding the catalase KatB) was induced at 18°C, which 
would be inconsistent with our hypothesis about an in- 
active status for the Fur protein at low temperatures. 
However, several reports have indicated an association 
between decreased environmental temperature and in- 
creased intracellular oxidative stress in bacteria. Thus, 
there is evidence that free radical production (super- 
oxide 0~ 2 , hydrogen peroxide H 2 0 2 , or hydroxyl radical 
HO") in bacterial cells is stimulated at low temperatures, 
apparently in an iron-independent manner. Therefore, 
the expression of oxidative stress adaptation genes, such 
as catalases, increase considerably [48,49]. A similar re- 
sponse may occur in our strain, justifying the observed 
induction in the catalase gene, as low temperature in- 
duces free radical production in cells, in turn increasing 
catalase production. The expression of the gene encod- 
ing catalase (KatB) was evaluated by RT-PCR analysis 
(Figure 3). Furthermore, it has been reported that iron- 
starvation inducible genes are also induced in response 
to oxidative stress in P. aeruginosa. This response ap- 
pears to be due to a transient loss of Fur repressor func- 
tion [50]. These observations are consistent with our 
data and support our hypothesis about the inactive sta- 
tus of the Fur protein at low temperatures. 

Additionally, within Cluster 6, we also found PSPPH_ 
1309, which encodes the cysteine desulfurase IscS, and 
PSPPH_1311, which encodes iron-sulfur cluster assem- 
bly protein IscA, both components of ISC (iron-sulfur 
cluster) system essential in the biogenesis of iron-sulfur 
(Fe-S) proteins in bacteria. It has been observed that 
some pathways involved in Fe-S cluster assembly operate 
under iron starvation and oxidative stress conditions 
[51,52], which agrees with the results obtained. 

On the other hand, several reports have indicated a 
correlation exists between the uptake- transport iron sys- 
tem and motility process and biofilm formation. Thus, 
iron deficiency stimulates twitching motility, a form of 
surface motility that is inconsistent with microcolonies 
and biofilm formation [53]. This is consistent with the 
results obtained in our experiments, because iron metab- 
olism genes and siderophores production are induced, 
simulating iron deficiency conditions, and motility pro- 
cesses appear to be favored, whereas biofilm or extracellu- 
lar polysaccharide formation is decreased (see data below). 

Hypothetical proteins and proteins with unknown 
function are induced at 18°C 

Among the differentially regulated genes induced at 
18°C, we found 15 genes that hypothetically encode 



conserved proteins (Cluster 7). Additionally, Cluster 8 
has genes that could not be grouped into any specific 
biological process but showed high transcript levels at 
18°C relative to 28°C. Within this cluster are genes enco- 
ding various transcriptional regulators, a gene that encodes 
an ATP-dependent helicase, DinG family (PSPPH_1406), 
and the PSPPH_4151 gene that encodes RNA polymerase 
sigma-54 factor RpoN whose expression was validated by 
RT-PCR assays (Figure 3). 

Low temperature represses alginate synthesis in ft 
syringae pv. phaseolicola NPS3121 

With regard to repressed genes identified at 18°C, we 
found a group of genes involved in alginate synthesis 
(Cluster 9). The genes enconding AlgX (PSPPH_1112), 
AlgG (PSPPH_1113), AlgE (PSPPH_1114), AlgK (PSPPH_ 
1115), and AlgD (PSPPH_1118), as well as the PSPPH_ 
1119 gene that encodes a hypothetical protein, were 
included in this cluster. Alginate is an extracellular poly- 
saccharide (EPS) produced by bacteria that is secreted into 
growth media and involved mainly in biofilm formation. 
Production of this co-polymer by P. syringae and P. 
aeruginosa has been previously reported [54,55]. Alginate 
production by P. syringae has been associated with in- 
creased epiphytic fitness, resistance to desiccation and 
toxic molecules, and the induction of water-soaked lesions 
on infected leaves. Studies have shown that alginate func- 
tions in the virulence of some P. syringae strains and facil- 
ities the colonization and/or dissemination in plants [55]. 
Although P. syringae pv. phaseolicola virulence is favored 
by low temperature, alginate production by this strain ap- 
pears to be repressed under these conditions. RT-PCR 
analyses confirmed the repression mediated by low tem- 
peratures of algD, the first gene of the alginate biosyn- 
thetic operon (Figure 3). The repression of alginate genes 
mediated by low temperature also has been observed in 
P. syringae pv. syringae, where the expression of algD, 
was induced at 28°C and significantly lower at 18°C [56]. 
To validate the microarrays results in P. syringae pv. 
phaseolicola NPS3121, the effect of temperature on EPS 
production (including alginate) was evaluated. Quantita- 
tive analyses showed that at 18°C the production of EPS is 
lower (76.65 ± 4.09 \ig) compared to when the bacterium 
is grown at 28°C (192.43 ±14.11 |ig). Thus, the results 
demonstrate that the low temperatures decrease EPS pro- 
duction by the bacterium. Alginate gene regulation is 
complex and varies between species. In P. aeruginosa, it 
has been reported that sigma factor- 54 (RpoN) represses 
algD expression by sigma factor antagonism [57]. A simi- 
lar phenomenon could be occurring in our strain, because 
the expression levels of the rpoN gene (PSPPH_4151) 
are consistent with the low expression of alginate genes. 
Furthermore, it has been reported that a coordinated 
expression exists between flagellum synthesis and EPS 
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production. In P. aeruginosa, the FleQ protein, a master 
regulator of flagella genes, represses the expression of 
genes involved in EPS synthesis, leading to planktonic 
cells. When this repression is released, the flagellum genes 
are repressed and EPS production is favored [58]. The al- 
ginate gene repression observed in our microarray, could 
also be due to repression exerted by FleQ protein, which 
is induced in our experiment, in a similar manner to what 
occurs in P. aeruginosa. Thus, the results of the microarray 
are consistent with the fact that EPS production (e.g., 
alginate) is decreased at low temperatures whereas expres- 
sion of motility genes is favored. 

Genes involved in plant-pathogen interactions are 
affected by low temperature 

Another group of genes repressed at 18°C comprise 
Cluster 10 and correspond to genes related to plant- 
pathogen interactions (Table 2). Four of these genes 
encode Type III effector proteins (T3EFs): HopABl, 
HopWl, HopDl, and AvrB2, but the other genes are 
involved in cell wall degrading enzyme synthesis, such as 
pectin lyase (PSPPH_3992) and polygalacturonase (PSPPH_ 
A0072). The repression of some T3EFs genes and cell 
wall degrading enzyme genes was validated by RT-PCR 
(Figure 3). The classification of these genes as pathogenicity 
and/or virulence factors in P. syringae pv. phaseolicola has 
been previously reported [18]. It known that phytopatho- 
genic bacteria suppress plant innate immunity and promote 
pathogenesis by injecting directly into host cells effector 
proteins (T3EFs) by a type III protein secretion system 
(T3SS). However, in the majority of cases, neither the mode 
of action nor the targets of these effector proteins within 
the plant are known [59]. In addition, phytopathogens 
synthesize and secrete various cell wall degrading enzymes, 
which facilitate pathogen entry and nutrient release for its 
growth [1]. Thermoregulation of these genes has been ob- 
served in other bacterial phytopathogens, where their ex- 
pression is favored at low temperatures, a phenomenon 
opposite to data obtained in our experiments [4]. However, 
it has been reported that in P. syringae pv. tomato DC3000, 
iron bioavailability regulates the expression of T3SS com- 
ponent genes. Thus, high iron concentrations induced ex- 
pression of genes such as hrpRS and hrpL, which in turn 
regulates the expression of T3SS genes, by an as yet 
unknown mechanism [60]. Based on this, our microarray 
results might be explained by the fact that the uptake- 
transport iron genes were induced, mimicking iron limiting 
conditions, which could lead to the observed repression of 
T3SS genes. 

Genes related to the Type IV secretion system (T4SS) are 
repressed at 18°C 

Another group of genes differentially repressed at 18°C 
comprise Cluster 11. They include genes related to the 



type IV secretion system (T4SS), which is closely related 
to systems involved in the conjugal transfer of DNA 
(Table 2). Nine of these genes encode conjugal transfer 
proteins and two encode transcriptional regulators, all 
within plasmid B (pPhl448B) of P. syringae pv. phaseolicola 
(Figure 2) [18]. The pPhl448B plasmid belongs to the well- 
described pPT23A plasmid family, whose members have 
been demonstrated to play an important role in the inter- 
action of the P. syringae pathogen with host plants [61]. 
Many of the pPT23A family plasmids are known to be 
conjugative plasmids. The putative T4SS encoded in the 
pPhl448B plasmid has been classified as a type IVA system, 
due to its high similarity with the type IV secretion genes of 
Agrobacterium tumefaciens (the virB operon and virD4) 
[61]. In addition to functioning in conjugation, some bac- 
terial T4SSs are capable of delivering effector proteins or 
toxins to host cells, thus acting as virulence factors [62]. 
However, in the case of P. syringae pv. phaseolicola 1448A 
T4SS, it has been suggested to have a role in conjugal trans- 
fer of DNA rather than virulence-related protein transloca- 
tion [18]. Thermoregulation of some T4SS genes in various 
bacteria has already been reported, similar to our results in 
this study [4]. More experimental work is necessary to elu- 
cidate the role of these genes in P. syringae pv. phaseolicola 
NPS3121 and their relationship to temperature. 

Low temperature represses the heat shock response 

Another group of genes repressed at 18°C correspond to 
those encoding heat-shock proteins (Cluster 12). Genes 
that encode the HslVU and GrpE heat-shock proteins, 
as well as the genes encoding DnaK, GroEL, and ClpB 
chaperones were included in this cluster. Heat shock 
proteins (HSPs) are a class of functionally related pro- 
teins that are responsible for monitoring the state of 
protein folding in cells. They function as molecular 
chaperones, facilitating the folding of partially or fully 
unfolded proteins. Their expression is increased when 
cells are exposed to elevated temperatures or other 
stresses, to cope with protein damage. If however, the 
temperature decreases, a reverse response is observed 
and heat-shock gene transcription decreases [63]. This 
latter behavior is similar to the results obtained in our 
experiments, where the low temperature decreased the 
transcript levels of heat-shock genes. In E. coli, HSP syn- 
thesis is repressed during growth at low temperatures 
[64]. A similar response has been observed in P. putida, 
where low temperatures also decrease the expression of 
these genes [65]. 

Transcription and replication are repressed by low 
temperature 

Cluster 13 includes genes involved in nucleic acid syn- 
thesis. Two of these genes (PSPPH_4598 and PSPPH_ 
4599) encode RNA polymerase beta subunits involved in 
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mRNA synthesis. Three of these genes (PSPPH_2495, 
PSPPH_B0043, and PSPPH_A0002) are related to the 
replicative process of DNA synthesis. This result suggests 
that both processes are affected by low temperature in 
P. syringae pv. phaseolicola NPS3121, which is consistent 
with the decreased growth rate observed. This behavior is 
similar what was previously observed in P. putida where 
low temperature also reduces proteins involved in the 
transcription and replication processes [65]. 

Finally, similar to the analysis and clustering of acti- 
vated genes, repressed genes at 18°C that hypothetically 
encode conserved proteins were grouped into Cluster 
14. Likewise, those genes whose products could not be 
grouped into any specific biological process were in- 
cluded in Cluster 15. The relationship of these genes to 
the physiology of the bacterium to low temperatures re- 
mains unknown and more experimental work is required. 

Conclusions 

In general, the results of the microarray provided us 
with a global view regarding the physiology of P. syringae 
pv. phaseolicola NPS3121 at 18°C, a temperature that 
favors disease development and produces the major viru- 
lence factor of the bacterium. The results showed that 
common processes in response to low temperature, such 
as cell-envelope remodeling, transcription, translation, 
and the heat-shock response, are also affected in this bac- 
terial phytopathogen. In addition, low temperatures influ- 
ence phaseolotoxin synthesis as well as the expression of 
various virulence factors involved in disease development. 
Furthermore, our data show low temperature-dependent 
expression of T6SS, thus being the first report about the 
expression of this cluster of genes in P. syringae pv. 
phaseolicola. In general, the expression profile obtained in 
this study suggest that low temperatures generate an oxi- 
dative stress in the bacterium, which leads to expression 
of uptake-transport iron genes (simulating iron starvation 
conditions) that in turn are related to the expression of 
various processes such as motility, biofilm production, 
and T3SS. From the data obtained in this study, we can 
begin to understand the temperature dependent strategies 
used by this phytopathogen during host interactions and 
disease development. 

Methods 

Bacterial growth conditions and RNA isolation 

The P. syringae pv. phaseolicola NPS3121 strain was 
grown at 18°C and 28°C in M9 minimal media sup- 
plemented with 0.8% glucose as the carbon source. 
The growth conditions were as follows: pre-inoculums 
(25 mL) of P. syringae pv. phaseolicola were grown in 
M9 minimal media overnight at 28°C. The cells were 
washed once with M9 media and inoculated into 200 
mL of M9 minimal media at optical density (OD 6 oonm) 



0.1. To evaluate the effect of temperature, the cultures 
were incubated at 18°C or 28°C and grown until they 
reached the transition phase (OD 600nm 1.1 at 18°C and 
1.2 at 28°C) and RNA was extracted. For RNA isolation, 
the cells were recovered by centrifugation at 10,000 rpm 
for 10 min at 4°C, washed with sterile deionized water, 
and stored at -80°C. The supernatants from each culture 
were removed for phaseolotoxin production assays. 
Total RNA was extracted using the TRIzol Reagent fol- 
lowing the manufacturer's instructions (Invitrogen, CA, 
USA). A second purification step was performed using 
RNeasy MinElute spin columns (Qiagen, CA, USA) 
to remove any contaminating DNA. RNA was eluted in 
50 uL of diethylpirocarbonate (DEPC) -treated water. 
RNA concentration was determined using a ND-1000 
spectrophotometer (NanoDrop). RNA integrity was veri- 
fied by analytical agarose gel electrophoresis. 

Phaseolotoxin assays 

Phaseolotoxin production by P. syringae pv. phaseolicola 
was assayed using the E. coli JM103 strain growth inhi- 
bition assay as previously described [66]. In each case, 
plates containing arginine (10 mM) were used as con- 
trols to confirm that growth inhibition was due to 
phaseolotoxin effects. 

Microarray processing, data acquisition, and statistical 
analyses 

Previously, our group constructed a DNA microarray of 
P. syringae pv. phaseolicola NPS3121 (GEO accession 
number GPL7115) with 1911 probes of an average length 
of 2.4 kbp, which represents approximately IX of the 
P. syringae pv. phaseolicola NPS3121 genome. This micro- 
array contains also several PCR products corresponding 
to various genes with known functions that were printed 
as controls [67]. To perform this study, we used this 
P. syringae pv. phaseolicola NPS3121 DNA microarray. 
Each microarray experiment was repeated six times: two 
technical replicates with the same RNA samples and three 
biological replicates using RNA isolated from a different 
culture. cDNA synthesis, labeling, hybridization, washing, 
and chip scanning were performed at the Microarray 
Core Facility at CINVESTAV-LANGEBIO. Hybridized 
microarray slides were scanned (GenePix 4000, Axon In- 
struments, Inc) at a 10- um resolution, adjusting the laser 
and gain parameters to obtain similar levels of fluores- 
cence intensity in both channels. The spot intensities were 
quantified using Axon GenePrix Pro 6.0 image analysis 
software. First, an automatic spot finding and quantifica- 
tion option of the software was used. Subsequently, all 
spots were inspected individually and in some cases, the 
spot diameters were corrected or the spots were removed 
from the analyses. The mean of the signals and the me- 
dian of backgrounds were used for further analyses. Raw 
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data were imported into the R.2.2.1 software. Background 
signals were subtracted using Robust Multichip Analysis 
(RMA) whereas the normalization of the signal intensities 
within slides was carried out using "print-tip loess" 
method and the LIMMA package. Normalized data were 
log2 transformed and fitted into mixed model ANOVAs 
using the mixed procedure. The ^-values of the low 
temperature (18°C) effect were adjusted using the False 
Discovery Rate method (FDR). Differentially expressed 
genes were identified using cut-off criteria of > 1.5 for up- 
regulated and <0.6 for down-regulated genes (FDR, 
p-value < 0.05). Analyses of distribution and the loca- 
tion of differentially expressed genes in the P. syringae 
pv. phaseolicola 1448A sequenced genome were 
performed using the GenoMap software [68]. 

Microarray validation by reverse transcription-PCR 
analyses 

To validate the microarray data, we performed reverse 
transcription (RT)-PCR analyses. The expression levels 
of several genes with different biological functions were 
evaluated by this technique. These experiments involved 
independent biological experiments from those used for 
microarray analyses. DNA-free RNA was obtained as de- 
scribed above and the integrity of the RNA was evalu- 
ated by agarose gel electrophoresis. Total RNA (200 ng) 
was used for RT-PCR using the Superscript one-step kit 
(Invitrogen). Controls used for each set of primers were 
1) PCR without the reverse transcription step to verify 
the absence of DNA, 2) RT-PCR performed without 
RNA templates to detect any contaminating DNA/RNA, 
3) PCR performed using genomic DNA as a template to 
ensure primer fidelity, and 4) amplification of a portion 
of the 16S rRNA operon as an internal control of the 
reaction. The RT reaction was performed at 50°C for 
30 min. PCR amplification was performed at 94°C for 
2 min for 1 cycle; 94°C for 30 s, 55-58°C for 30 s, and 
72°C for 1.0 min for 20-28 cycles; and 72°C for 10 min 
for 1 cycle . 

Molecular biology techniques 

Routine techniques were performed using standard 
protocols [69]. Genomic DNA of P. syringae pv. 
phaseolicola NPS3121 was isolated as described previously 
[70]. PCR products were amplified with Platinum supermix 
(Invitrogen). Primers were designed using Vector NTI 
Software (Invitrogen), with reference to the previously 
reported sequence of the 1448A strain (Gene Bank acces- 
sion no. CP000058) [18]. The oligonucleotide primers used 
in this study are listed in Additional file 1. 

Motility assays 

To evaluate the motility of P. syringae pv. phaseolicola 
NPS3121 and the influence of temperature on this process, 



three strategies were used. The swimming and swarming 
motility of P. syringae pv. phaseolicola NPS3121 were 
assessed on semisolid KB plates containing 0.3% and 0.5% 
agar, respectively, as described in previous studies [41,42]. 
The cells were grown in KB broth overnight at 28°C, and 
harvested and resuspended in KB to OD 600 = 1. 50 uL 
of bacterial suspensions were inoculated on filter disks 
(6 mm in diameter) and placed in the center of the plate. 
Plates were incubated for 24 h at 28°C and 18°C before 
photography. A second strategy was performed to evaluate 
the swimming and swarming motility of P. syringae pv. 
phaseolicola NPS3121. To ensure that the bacteria were in 
the same physiological condition as when the transcrip- 
tome analysis was performed, the P. syringae pv. 
phaseolicola NPS3121 strain was grown in M9 media at 
28°C and 18°C until they reached the transition phase. 
Bacterial suspensions (50 uL) were inoculated on filter 
disks (6 mm in diameter) and placed in the center of semi- 
solid M9 plates containing 0.3%, 0.4%, and 0.5% agar. 
Plates were incubated for 48 h at 28°C and 18°C. Finally, 
motility was also evaluated using the stab technique in 
semisolid KB and M9 media (0.3% and 0.5% agar) in 
glass tubes. The tubes were incubated at 28°C and 18°C 
for 48 h. As controls, we used the P. syringae strains 
pv. tomato DC3000 and pv. tabaci PTBR2004. Experi- 
ments were performed three times with three replicates 
per treatment. 

Quantification of siderophores 

Siderophore production into the culture supernatant by 
bacterial strains was determined using chrome azurol S 
(CAS) liquid assays as previously described [71]. Briefly, 
the P. syringae pv. phaseolicola NPS3121 strain was 
grown in M9 media at 28°C and 18°C until they reached 
the transition phase. The supernatant was recovered by 
centrifugation at 8,000 rpm for 15 min at 4°C and fil- 
tered through a 0.45-um-pore-size filter (Millipore). For 
siderophore quantification, a standard curve was prepared 
with desferoxamine mesylate. Experiments were performed 
three times with four replicates per treatment. 

Quantitative analyses of EPS 

The quantification of EPS was performed as previously 
described with some modifications [72,73]. The P. 
syringae pv. phaseolicola NPS3121 strain was grown in 
M9 media at 28°C and 18°C until they reached the tran- 
sition phase [the growth stage in which the microarrays 
analysis was performed and the repression of EPS syn- 
thesis genes (alginate) was observed]. The bacterial 
cells were harvested by centrifugation at 8,000 rpm for 
15 min at 4°C. After centrifugation, the supernatant 
was mixed with three volumes of ice-cold 95% ethanol 
(with stirring) for 24 h at -20°C to precipitate the ex- 
tracellular polysaccharide (EPS). EPS was recovered by 
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centrifugation at 10,000 rpm for 20 min at 4°C. The pel- 
let was washed twice with 95% ethanol and once with ab- 
solute ethanol. Quantification of the EPS was performed 
using the phenol-sulfate method. Total EPS was measured 
using a glucose standard curve. Experiments were performed 
three times with four replicates per treatment. 
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Additional file 1: This Word file contains the sequence of 
oligonucleotides used in the RT-PCR assays. 
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